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Abstract 

The hydrodynamic force on a cylinder oscillating transversely to its axis is a nonlin- 
ear function of the displacement amplitude xo. We report measurements and numerical 
calculations of the force at frequencies low enough that 6 > R, where 6 is the viscous 
penetration length and R is the cylinder radius. For small amplitudes, the numerically 
calculated Fourier transform of the force per unit length, F sma ii, agrees with Stokes’ 
analytical calculation. For larger amplitudes, the force per unit length found by both 
calculation and measurement is F = F sraa uC {xq/6, R/6). The complex function C 
depends only weakly on R/6, indicating that xq/6 is more appropriate as a scaling 
variable than the Keulegan-Carpenter number KC = ttxq/R. The measurements used 
a torsion oscillator driven at frequencies from 1 to 12 Hz while immersed in dense 
xenon. The oscillator comprised cylinders with an effective radius of R = 13.4 fi m 
and oscillation amplitudes as large as xq/6 = 4 (corresponding to KC as large as 71). 
The calculations used similar conditions except that the amplitudes were as large as 
xq/6 = 28. 


NASA/CR— 2007-2 1 5050 


1 



1 Introduction 


Small transverse oscillations of a cylinder will generate a flow that is periodic and effectively 
two-dimensional. Describing the resulting force on the cylinder requires two independent 
variables because the oscillation frequency / = tv/2n can be varied independently of the 
velocity. The most frequent choices for the two variables are the Keulegan-Carpenter number, 

KC = tt|, (1) 

which compares the oscillation amplitude xq to the cylinder radius R , and the Stokes viscous 
parameter, 


0 


4 

7 r 



( 2 ) 


which compares the cylinder radius to the viscous penetration length defined by 


6 = 



( 3 ) 


where p and p are the viscosity and density. 

The chief result of this paper is that, when 6 > R (or (3 < 1), the scaling variable most 
appropriate for describing the hydrodynamic nonlinearity is the reduced length 

( 4 ) 
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This choice minimizes the dependence on the second variable. We support this result with 
a simple scaling argument, numerical calculations, and experimental measurements. Figure 
1 indicates the previously unexplored region of KC and 3 spanned by the present measure- 
ments and shows the locations of several familiar oscillators. 

This work was made in preparation for an experiment, called “Critical Viscosity of 
Xenon” (CVX-2), that measured the shear thinning predicted to occur near the liquid- 
vapor critical point [15]. Shear thinning caused the response of the CVX-2 viscometer to 
be nonlinear, but the convective term of the Navier-Stokes equation caused a much larger 
nonlinearity [16]. To understand the larger nonlinearity we calculated and measured the os- 
cillator’s response when the xenon was far from the critical point and had no shear thinning. 

CVX-2 flew on board Space Shuttle Columbia in 2003. The microgravity in Earth orbit 
was necessary to obtain the homogeneity of density and temperature sufficient to observe 
shear thinning. For the forseeable future, the CVX-2 results will likely be the only data on 
shear thinning in a simple pure fluid. Correctly interpreting those data, as well as addressing 
the classic problem of an oscillating cylinder, motivated the present study. 

The experimental frequencies used by CVX-2 (1 to 12 Hz) corresponded to 6 > R (or 3 < 
1). With one exception [11], previous investigations dealing with the force on the cylinder, 
the surrounding flow patterns, or both [1, 2, 3, 4, 5, 6, 7, 8, 9, 10] have examined flows with 
KC similar to or smaller than those used here, but the values of 3 were always much larger 
than 1. (The steady streaming induced by an oscillating cylinder is not addressed here.) 
Large values of 3 occur in practical applications such as pilings for oil drilling platforms, and 
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they are compatible with the usual boundary-layer assumption that b <y; R. In contrast, the 
present results are in the range 0.01 < f3 < 0.2, which means that 6 > R. 

The exception is the numerical work by Morfey and Tan [11], which examined the range 
10 G < P < 0.3. Their calculations were motivated by the attenuation of large-amplitude 
sound in fibrous blankets. More discussion? 

Figure 1 can be related to the flow visualization study by Tatsuno and Bearman [10]. 
In the range defined by 1.6 < KC < 15 and 5 < (3 < 160 they found eight regimes of flow 
patterns. Only two regimes, which they called A and A *, were two-dimensional. Figure 2 
shows that the curve xq/6 = 7 approximately bounds regimes A and A*. As shown in Figure 
1, an extrapolation of that boundary to small 3 indicates that the present experimental 
results were obtained within the two-dimensional regime. 

This paper studies the viscous force per unit length F (t) that opposes the motion of a 
cylinder whose time dependent displacement is 

x ( t ) = Re [x 0 e ja;t ] . (5) 

For small amplitudes, the Navier-Stokes equation reduces to the unsteady Stokes equation, 
and the force is periodic and proportional to xq. The proportionality allows the F (t) to be 
written as the product of the fluid mass per unit length displaced by the oscillator n R 2 p, 
the oscillator’s acceleration amplitude cu 2 x o, and a dimensionless reduced force B (R/6) as 


NASA/CR— 2007-2 1 5050 


4 



follows [12]. 


F (t) = Re 


small 


(/) e 


iujt 


= Re 


-i (irR 2 p(jJ 2 x 0 ) B 


Atot 


( 6 ) 


This form is convenient for an oscillator used as a viscometer because 8 can be varied by 
changing the frequency while the oscillator is immersed in a fluid of known density and vis- 
cosity. The choice of R used in B ( R/8 ) is arbitrary; therefore, measuring the Fourier trans- 
form . F STr i a .ii (/) while varying the frequency yields the calibration function B (R/8) without 
knowledge of the oscillator’s geometry. Figure 3 shows the magnitude and phase of B (R/8) 
measured for the present experimental oscillator at small amplitudes [13]. Also shown is 
Stokes’ solution for a cylinder with a circular cross section. Its modern formulation [14] is 


Be yl 



(k' + ik) = i 


\ 4 KM 

. zK 0 (z)_ 


( 7 ) 


where K n (z) is the modified Bessel function of order n, and z = (l + i) (R/8). Here and in 
the remainder of this paper, we will use R/8 instead of ft to make clear the role of 8. 

For moderate amplitudes, the flow remains periodic, but the force on the cylinder is no 
longer proportional to Xq. Therefore we generalized Eq.(6) by multiplying the reduced force 
B by a complex function C that depends on reduced amplitude Xq/8 as well as R/8. 




( 8 ) 
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The amplitude function C, which describes our numerical and experimental results, is written 
as a function of x 0 /6 instead of KC . Figure 4 illustrates how choosing x$/8 for the first 
variable minimizes the dependence on the second variable R/6. The magnitude and phase 
of C numerically calculated for a circular cylinder are each described approximately by a 
single function of Xq/8. 

The following section describes a simple analytical argument that supports the impor- 
tance of Xq/8 as a scaling variable. Following that is a section that describes numerical 
calculations for the force on an oscillating circular cylinder. The final section compares the 
experimental measurements to the numerical results. 

2 Scaling of the Navier- Stokes equation 

Imagine a solid body that oscillates in an incompressible fluid according to x (t) = (Co cos (cut). 
For small amplitudes xq, the velocity u and pressure p will be periodic everywhere in the 
fluid. For moderately larger Xq, the velocity and pressure will remain periodic, even though 
harmonics and steady streaming may appear. Because there is a obvious time scale a;" 1 one 
can define dimensionless time, velocity, pressure, and spatial derivatives as follows. 


t. M = cut 

(9a) 

u* = u/ (uuxq) 

(9b) 

L 

(9c) 

^3 

* 

III 

^3 

rjcuxo 

V* = LV and = L 2 V 2 

(9d) 
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Here * indicates the dimensionless quantities. The characteristic length L is not yet iden- 


tified, which precludes its use in a viscous time scale such as pL 2 /^; doing so leads to an 
ambiguous result. In terms of these dimensionless quantities, the Navier-Stokes equation is 

<9u* x Q , _ . 8 _ 8 _ 

^- + — (u*- V*)u* = ^— V*p* + — V*u*. (10) 

This equation can be expressed in terms of a single dimensionless scaling variable by requiring 
that 


L = Sf L (^) , (11) 

where /l{xq/8) is an arbitrary function of x 0 /8. The simplest choice, L = 8, puts the 
Navier-Stokes equation into the form 

^ + (y ) (u* • V,) u* = 1 [- + VX] • (12) 

Eq.(12) shows that the Navier-Stokes equation for a periodic flow can be scaled by the 
reduced amplitude Xq/8 if the scales chosen for time and length are respectively to and 
8. Note that Xq/8 is equivalent to an “oscillatory Reynolds number” whose characteristic 
velocity and length are respectively ojxq and 8/2. Similar Reynolds number have been used 
in discussions of streaming from oscillating bodies [17, 18]. 
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3 Calculations 


3.1 Numerical model 

The numerical model calculated the force on the experimental oscillator, which was con- 
structed from a screen comprising many wires of noncircular cross section. (See Figure 5 
and the Measurements section.) Four simplifications were used. First, we assumed that 
the oscillator behaved as an assembly of isolated wires and therefore could be represented 
by a single cylinder oscillating transversely to its axis. This approximation was reasonable 
because the spacing between adjacent wires was much larger than the viscous penetration 
length. Second, we assumed that the cross section of the wire could be represented by a 
circle of appropriate diameter. Although Figure 5 shows that the actual cross section of the 
screen wire was closer to a rectangle, numerical studies, illustrated in Figure 6, showed that 
the effect of the cross-section shape was small. These first two simplifications reduced the 
model to the classical flow problem of a rigid circular cylinder, surrounded by an otherwise 
stationary incompressible Newtonian fluid, and oscillating sinusoidally in a direction normal 
to its axis. See Figure 7. 

The third simplification assumed that the flow held was two-dimensional. An example of 
three-dimensional how is Tatsnno and Bearman’s regime B , where counter- rotating vortices 
appear alternately along the cylinder axis [10]. The fourth simplification assumed that 
the how held had reflection symmetry about the oscillation direction. Both simplifications 
were consistent with the patterns seen in Tatsuno and Bearman’s how regimes A and A*, 
although the extrapolated boundary in Figure 1 suggests that real hows with Xq/8 > 7 may 
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be three-dimensional. 


The values of density (p = 1116 kg-m~ 3 ) and viscosity (77 = 52 pPa-s) were similar 
to those of the experimental measurements. (Some test calculations used 7 / = 72 pPa-s.) 
The value of the cylinder radius ( R = 13.4 pm) was similar to the half width of a single 
screen wire. The range of oscillation frequencies, 0.5 < / < 12 Hz, included the experimental 
frequencies, but the largest amplitude of xq = 1000 pm exceeded the experimental amplitude 
at all frequencies. 

We used the moving reference frame approach described in Appendix A to handle the 
oscillatory motion of the cylinder. This approach transforms the original moving boundary 
problem into a fixed domain problem in which the cylinder is stationary and surrounded 
by oscillating fluid. To find the force F ( t ) on the cylinder in the laboratory reference 
frame, we first computed the corresponding force F' (t) in the moving reference frame. We 
then removed the spurious force due to the non-inertial frame by applying the correction 
prescribed in Eq.(33). The moving reference frame eliminates the need for re- meshing at 
each time step. In a stationary reference frame, the re-meshing cannot prevent distortion of 
the mesh, and the distortion at large oscillation amplitude makes convergence of the solution 
difficult. 

3.2 Numerical Techniques 

The fourth simplification mentioned above assumed that the solution had reflection symme- 
try about the laboratory x-axis passing through a diameter of the cylinder. Consequently, 
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the computational flow domain included only the upper half of the actual flow field defined 
by the width — w < x' < w, the height 0 < y' < h, and the top half surface of the circle. 
The computational domain and the boundary conditions are depicted in Figure 7. The size 
of the computational domain was w = h = 2 000 ft for all of the oscillation frequencies and 
amplitudes. 

The governing equations in the moving reference frame were solved using POLYFLOW 
[19], a commercial finite element method (FEM) program designed primarily for laminar, vis- 
cous flows involving complex rheology and free surfaces. Details of the FEM formulation and 
numerical techniques used in POLYFLOW are documented in [20]. Galerkin’s method was 
adopted in the FEM discretization for the momentum equations, and the two-dimensional 
FEM mesh was built with 9-node, quadratic, quadrilateral elements, in which the veloc- 
ity was approximated by quadratic shape functions. The pressure was approximated as 
piecewise-linear and discontinuous on inter-element boundaries. 

No-slip boundary conditions were used at the cylinder surface. Fine grading of the mesh 
was necessary near the cylinder due to the strong radial dependences of velocity and stress in 
that region. Numerical studies of those dependences led us to give the FEM mesh elements 
next to the cylinder the small size Ai = R/ 158. The mesh grading away from the cylinder 
surface was controlled by fixing the ratio q of successive element sizes, defined as 

q = Aj/A j+ i = 0.6. (13) 

Here A j is the characteristic length of element j . and j increases with radial distance from 
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the cylinder. Figure 8 illustrates the mesh grading. 

The transient problem was solved by a predictor-corrector time integration scheme in 
which the second-order Crank-Nicholson method was selected for the corrector. At each time 
step, the non-linear algebraic system resulting from the FEM discretization was solved by 
the Newton-Raphson iteration scheme. Termination of the nonlinear iteration was controlled 
by an iteration convergence tolerance of 10 4 for the relative error norms of the residuals of 
the governing equations. For each simulation, we calculated the transient solution for 600 
time steps, yielding a total of 6 complete oscillation cycles. Steady state was achieved after 
4 complete oscillation cycles, so the fifth and sixth oscillation cycles were used to determine 
the Fourier coefficients of the force on the cylinder. 

The velocity u 0 = xocu sin (cut) imposed at the domain boundaries was sinusoidal at 
frequency cu. However, the resulting force on the cylinder F (t), after correcting for the 
moving reference frame, contained harmonics for large amplitudes Xq. 

OO 

F (■ t ) = a 0 + [ b n cos (nut) — a n sin (ncut )\ , (14) 

n= 1 

Several recent analyses [21, 22, 23, 24, 25, 26] of large-amplitude shear viscometry decom- 
posed the non-linear force into a series of harmonic functions in a similar manner. The 
discrete Fourier analysis of F ( t ) was transformed into a general linear least-squares curve 
fitting problem and solved using the singular value decomposition algorithm described in 
[27]. A total of 41 terms in Eq.(14) were used to fit the Fourier coefficients a n and b n with 
n = 0, 1, • • • , 20. Numerical tests indicated that the accuracy of the Fourier coefficients 


NASA/CR— 2007-2 1 5050 


11 



fitted by the singular value decomposition algorithm was satisfactory. This paper discusses 
the results at only the fundamental frequency (n = 1) because the harmonics were small. 
The magnitude at 3a; never exceeded 7 % of that at cu, and the magnitudes of higher odd 
harmonics were 1 % or smaller. The magnitudes for all even harmonics (2cu,4cu, ...) were 
insignificant. 

3.3 Validation 

When £0 <C R, the force on the oscillating circular cylinder can be calculated with Stokes’ 
analytical solution. If the cylinder moves according to x (t) = Xq cos (cut), the opposing force 
per unit length is 


F (t) = — a\ sin (cut) + b\ cos (cut) = Re (b\ + ia\) e iuJt J , (15) 


where 


rf^ =Re [^‘] “ d ^JpAo = ~ lm[B ^' V (16) 

We estimated the accuracy of the numerical solutions at small amplitude by comparing 
them with the analytical solutions. Table 1 shows example results for three values of the 
mesh grading ratio q. The value q = 0.60 reduced the errors for a± and b\ to less than 0.06 % 
of the magnitude for all frequencies. (Fortuitous cancellations made the magnitude errors 
even smaller.) 
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Table 1. Small-amplitude comparisons between Stokes’ analytical solution and 
numerical solutions using three mesh grading ratios q. The Fourier coefficients 
Oi and bi of the force on the circular cylinder were calculated for p = 1116 kg/m 3 
and p = 72 /iPa-s at / = 1 Hz. The cylinder radius and oscillation amplitude 
were 13.4 pin and 1 pm respectively. The errors are relative to the magnitude 

\A*i + b\. 


mesh 


error 


error 


error 

q = 0.80 

-2.38163 

1.18 % 

0.835909 

1.01 % 

2.52407 

1.45 % 

q = 0.65 

-2.35201 

-0.01 % 

0.810190 

-0.03 % 

2.48764 

-0.02 % 

q = 0.60 

-2.35270 

0.02 % 

0.810376 

-0.02 % 

2.48836 

0.01 % 

Stokes solution 

-2.35227 


0.810815 


2.48809 



The effect of mesh grading was studied at large amplitude also. Table 2 compares the 
result for q = 0.60 to those obtained for coarser meshes. As with the small-amplitude 
comparisons, the results are nearly independent of q for q < 0.65. The present results were 
obtained with q = 0.60, so the uncertainty due to that mesh size is likely less than 0.1 %. 

Table 2. Large-amplitude comparisons between numerical solutions using three mesh 
grading ratios q. The Fourier coefficients (i\ and b\ of the force on the circular cylinder were 
calculated for 77 = 52 //Pa-s at / = 12 Hz and xq = 1000 pm. The differences from the 
q = 0.60 result are relative to the magnitude \j + b\ . 


mesh 

aiXlO 6 

difference 

biXlO 6 

difference 

y/af + xlO 6 

difference 

q = 0.80 

-139.0283 

8.76 % 

43.01326 

8.92 % 

145.5301 

10.82 % 

q = 0.65 

-127.4870 

-0.03 % 

31.14365 

-0.13 % 

131.2359 

-0.06 % 

q = 0.625 

-127.4856 

-0.03 % 

31.21592 

-0.07 % 

131.2517 

0.05 % 

q = 0.60 

127.5294 


31.30997 


131.3166 



The effect of domain size was studied by truncating the domain lengths L and H from 
2000.R to 500 R. At / = 1 Hz and x 0 = 1000 pm, the reduction changed cp and bi by only 


0.02 %. 
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3.4 Numerical results 


Figure 9 shows examples of the flow held near the cylinder at three oscillation amplitudes 
Xq. The contours show constant values of the shear component of the rate-of-deformation 
tensor, defined as D xy = ( du/dy ' + dv/dx ') /2, where u and v are the velocity components 
in the x' and y' directions. 1 

Tables 3 and 4 give the magnitude and phase of the Fourier-transformed force per unit 
length F (/) = l>i + ia± normalized by 7rR 2 puJ 2 x 0 . The ranges of frequency and amplitude 
include the experimental values. 

Table 3. Magnitude of F (/) / {jiR 2 poj 2 x 0 ) computed for R = 13.4 pm, p = 1116 
kg-rn 3 , p = 52 pPa-s. 


x 0 (pm) 

0.5 Hz 

1 Hz 

3 Hz 

5 Hz 

8 Hz 

12 Hz 

1 

134.499 

77.407 

33.529 

23.180 

16.725 

12.766 

10 

134.511 

77.422 

33.549 

23.203 

16.751 

12.795 

60 

134.942 

77.940 

34.233 

23.972 

17.596 

13.700 

100 

135.714 

78.849 

35.360 

25.174 

18.836 

14.942 

200 

139.056 

82.545 

39.199 

28.850 

22.267 

18.121 

400 

149.271 

92.251 

47.068 

35.770 

28.386 

23.638 

600 

160.662 

101.913 

54.122 

41.842 

33.702 

28.405 

800 

171.754 

110.942 

60.534 

47.337 

38.501 

32.699 

1000 

182.308 

119.391 

66.477 

52.429 

42.955 

36.692 
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Table 4. Phase of F (/) in radians. 


x 0 (Aim) 

0.5 Hz 

1 Hz 

3 Hz 

5 Hz 

8 Hz 

12 Hz 

1 

-1.2592 

-1.2196 

-1.1407 

-1.0964 

-1.0509 

-1.0079 

10 

-1.2593 

-1.2198 

-1.1414 

-1.0976 

-1.0530 

-1.0111 

60 

-1.2626 

-1.2268 

-1.1638 

-1.1346 

-1.1093 

-1.0889 

100 

-1.2682 

-1.2383 

-1.1945 

-1.1778 

-1.1635 

-1.1496 

200 

-1.2901 

-1.2768 

-1.2623 

-1.2520 

-1.2365 

-1.2166 

400 

-1.3391 

-1.3387 

-1.3249 

-1.3078 

-1.2846 

-1.2596 

600 

-1.3736 

-1.3726 

-1.3526 

-1.3329 

-1.3097 

-1.2880 

800 

-1.3963 

-1.3934 

-1.3701 

-1.3507 

-1.3298 

-1.3113 

1000 

-1.4122 

-1.4078 

-1.3836 

-1.3656 

-1.3468 

-1.3300 


We now use the amplitude function C (xq/6, R/6) to plot the numerical results. Com- 
bining Eqs.(6) and (8) gives 


C ^ |C| e* = 'j) , = (tl + iai) " „ 1 ^ . (17) 

lim x . 0 ^ 0 ( f / x o ) x ° Ph + ^Di^m 

Dividing the ratio F/x o calculated at amplitude .x (J by the ratio calculated at the smallest 
amplitude of 1 //m gives the magnitude and phase of C. Figure 4 plots | C\ and (j) as functions 
of reduced amplitude. The points have some scatter because R/6 varied by a factor of \[ 7 2A 
when the frequency was varied from 0.5 Hz to 12 Hz. The relatively small scatter shows that 
the dependence of C on R/6 is weak. 

We obtained a correction for the dependence of C on R/6 that allowed better comparisons 
between the results at different frequencies. This was done by reducing R from 13.4 /jm to 
6.7 /iin and repeating the 5 Hz calculations. The reduction of R caused small decreases of 
both \C\ and q i>, which allowed us to estimate the derivatives d \C\ / dR and d(j)/dR. Figure 
10 shows that plotting these derivatives as functions of x 0 /6 suggests empirical forms for the 
dependence of C on R/6. 
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The plot of d \C\ /dR showed that the dependence of \C\ on R/6 could be approximately 


described by 


1 d\C\ |C7| — 1 
| C\d(R/6) ~ 2 


Solving this differential equation yields 


R\ 

\8'6j 


1 - C 0 (Jj) exp 



-l 


(18) 


(19) 


where Cq(xq/8 ) is an undetermined function. Similarly, the plot of (j) showed that the 
dependence of (j) on R/6 could be approximately described by 


d(j) , 
d (R/6) ~ ^ 

Solving for (j) yields 

where (xq/8) is an undetermined function. 

Eqs.(19) and (21) use known functions of R/6 to describe the dependence of C on R/6. 
Therefore, the value of C calculated (or measured) for an oscillating cylinder of reduced 
amplitude xq/6 and reduced radius R/6 can be adjusted to the value of C expected at the 
same Xq/8 but different R/6. This adjustment is useful because it allows a calculation of C 


( 21 ) 
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obtained at one frequency to be compared to a calculation obtained at a different frequency 
by adjusting the calculations to a common “reference” value (R/8) re{ . The most convenient 
reference is (R/6) ief = 0, which corresponds to the limit of infinitesimally small radius or, 
equivalently, infinitesimally small frequency. The following equations adjust C (xq/8,R/8) 
to C (xq/8, 0). 


C 





1 

\C(x 0 /8,R/8)\ 

exp H) 


exp 



(22a) 

(22b) 


Figure 11 uses Eq.(22) to plot the magnitude and phase of C (xq/8,0) as a function of 
xq/8. Using (R/8) re{ = 0 improved the consistency of the results for xq/8 < 10, especially 
in the phase. For larger Xq/8 , the magnitude results are inconsistent by as much as 3 %. 
The cause of this remaining inconsistency may be due to the approximations in Eqs.(18) and 
(20). It is worth noting that, as indicated in Figure 1, the data for Xq/8 > 7 he outside the 
boundary extrapolated from Tatsuno & Bearman’s regimes A and A*. 

Note the initial quadratic behavior of \C\ and (j) in Figures 4 and 11. Such behavior is 
expected because an expansion of C about x 0 /8 = 0 must include only terms even in x 0 /8 ; 
otherwise the value of C would depend on the sign of Xq. The initial quadratic behavior 
was built into the curves drawn in Figure 11. These empirical descriptions of the numerical 
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results are given by Eq.(23), where b 0 = 4.3, c 0 = 36, a $ = 4, = 3.5, and = 36. 


C 




Pp c 0 

_(x 0 /S) 1/2 (xq/6) 2 _ 

n X- _i_ C <P 

0 (x 0 /<5) 1/2 0V<5) 2 


(23a) 

(23b) 


4 Measurements 


4.1 Experimental method 

The oscillator and xenon sample were the same as for the CVX experiment [13]. The density 
was fixed at xenon’s critical density (p = 1116 kg-m~ 3 ) by sealing the sample cell, and the 
viscosity was fixed near // = 52 //Pa-s by controlling the temperature near 295 K. A brief 
description of the apparatus is given below, and a detailed description can be found in [13]. 
A later publication will describe the electronic modifications made for CVX-2 that allowed 
operation at large Xq [16]. 

Figure 5 indicates the geometry of the oscillator. The oscillator was an 8 x 19 mm 
rectangle of screen that was cut out of a larger piece of nickel screen. The screen consisted of 
a square grid of wires formed by electrodeposition. When the screen was cut, one wire was 
left extending from the middle of both of the rectangle’s long edges. These two extensions 
served as torsion springs, and they were soldered to a stiff yoke that was centered between 
four electrodes placed parallel to the screen. While immersed in the xenon, the screen was 
driven to oscillate much like a child’s seesaw. 
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The viscous penetration length was always in the range 35 pm < 8 < 122 pm. Thus <5 was 
larger than the width of a screen wire (30 pm), but much less than the distance between wires 
(847 pm). Figure 5 shows that the cross section of a typical screen wire was approximately 
an 8 x 30 pm rectangle. Despite the wire’s noncircular cross section, Figure 3 shows that 
the drag on the screen resembled that of an assembly of isolated circular cylinders. Using 
the value R = 13.4 pm in Stokes’ analytical function B cy \ {R/8) approximately describes the 
experimental values obtained for B {R/8) at small amplitude. 

The oscillator was driven electrostatically at 1, 3, 5, 8, and 12 Hz. Four electrodes 
surrounding the oscillator (not shown in Figure 5) simultaneously applied combinations of 
a DC voltage Vbc and a sinusoidal voltage Vj n (t) = Vo sin {cut). The resulting torque on the 
oscillator was 


JVC*) 


Vdc 


■Hip {t) 

Xgap 


(24) 


where x t i p (t) is the displacement at the oscillator’s tip, and x gap ~ 4 mm is the length of the 
gap between the oscillator and one electrode. The dimensionless calibration factor K was 
adjusted so that the values of B {R/8) obtained for small x 0 agreed with those determined 
earlier for the CVX experiment [13], which also were made at small xq. 

The apparatus used for the present measurements differed from the CVX experiment in 
that the drive voltages Vbc an d VJ n it) were larger and the time dependence of Vj n it) was 
a single frequency instead of a frequency chirp. In particular, the second term of Eq.(24), 
which caused the electrostatic torque to increase linearly with the oscillator’s tip displacement 
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Xtip (t), was more significant. The result was a negative contribution to the oscillator’s spring 
constant, which the analysis handled by decreasing the oscillator’s resonance frequency by 
approximately 1 Hz from its vacuum value of 11 Hz. Terms of order (x t i p (t) /x gap ) 2 and higher 
were not included in Eq.(24) because the oscillator’s small amplitude and the approximate 
symmetry of the electrode pairs made them negligible. 

The oscillator’s time-dependent angular displacement 9 ( t ) was detected by the unbalance 
of a capacitance bridge operating at 10 kHz. The out-of-balance signal was fed to a lock-in 
amplifier, which generated a voltage whose oscillating component Vh ut ( t ) was proportional 
to 9 (t): 


9 (t) = 9 o cos (tut) 


1 

n i P 


f diCtip \ 

{dV^tJ 


Kut (t) • 


(25) 


Here, r t i p is the distance from the torsion axis to the oscillator’s tip. Video images of 
the oscillator were used to measure the derivative dx tip /dV out = (53.1 ± 1.2) /irn/V. The 
resolution of the images dominated the uncertainty of this result. 

The oscillator’s response was defined as the Fourier transform 9 (/) of the angular dis- 
placement divided by the Fourier transform N (/) of the applied torque. 


n ( f a \ = ill = ^ ( d Xti P \ ^out (/) 
1 ° j “ N(f)~r tip {dV ou J V„(/) 


(26) 


Here V- m (/) and V) U t (/) are the Fourier transforms of the input and output voltages. (For 
clarity, Eqs. (24) and (25) do not show the frequency-dependent corrections that account 
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for the anelasticity of the oscillator’s spring constant and the gain and phase shifts of the 
drive and detection circuits.) 

The amplitude Vo and frequency / of the AC drive voltage determined the oscillation 
amplitude Xq. A determination of the hydrodynamic force at amplitude Xq comprised two 
measurements. The first measurement used a drive torque sufficiently small that (7=1. 
(Its amplitude of 1.1 x ICC 10 N-m was approximately that used in the small-amplitude CVX 
experiment [13].) The second measurement used a torque that was larger by a factor of 
either 20 or 33 to produce the oscillation amplitude x 0 . Due to the oscillator’s response, the 
largest value of xq/S was achieved at 3 Hz. Figure 12 indicates the sizes of ay and b at 3 Hz 
relative to the wire cross section. 

The present measurements were made sufficiently far above the critical temperature of 
290 K that density stratification and internal waves were negligible. The viscosity was always 
within 1 % of the value used in the numerical calculations. 

4.2 Experimental analysis and results 

Each of the many wires comprised by the screen moved at a speed that depended on its 
distance from the torsion axis. We focused on the wire located at the oscillator tip, whose 
displacement was 


X t ip (t) = OqT tip COS (cut) = X 0 COS (cut) . 


(27) 


NASA/CR— 2007-215050 


21 



Comparing the experimental data to the numerical calculations at amplitude x (J required 
five steps. 

1. Eq.(34) in Appendix B was solved to convert the measured response G (/, 9 o) to values 
of B osc (R/6,9 q), which describes the hydrodynamics of the entire oscillator. Nonlin- 
earity of the capacitive drive and detection required a correction that did not exceed 
1.3 % of |G|. 

2. Eq.(37) in Appendix B was used to convert B osc ( R/6 , 9 0 ) and its derivative with respect 
to 9 0 into the product [BC\ that described the hydrodynamics of the oscillator tip. 
The value of B osc could not be compared directly to the numerical calculations because 
the total torque on the oscillator was the weighted sum of torques from all parts of the 
oscillator. 

3. The value of C at the oscillator tip was obtained by dividing [BC ] exp by B (R/6) as 
follows. 

r ,(x„R\ [B(R/S)C(x 0 /6,R/S) ]„ p 

wm (28) 

This approximated B (R/6) for a single wire by the experimental value for the entire 
oscillator B osc (R/6,9 q) measured at sufficiently small 9q. This approximation was 
accurate because the value of 9 0 used to determine B (R/6) corresponded to xq/6 < 0.1. 

4. The value of C was corrected for the noncircular shape of the wire’s cross section, 
which was modeled as a rectangle with quarter-circle ends (30 /iin width, 8 /iin radii). 
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The inset of Figure 6 indicates this shape, and the data show that the numerically 
calculated differences in C due to shape are small (about 1 % at the experimental 
maximum xq/8 = 4). The correction comprised interpolating functions fitted to the 
magnitude and phase data shown in Figure 6. The resulting values of C corresponded 
to the numerical calculations for a circle of 13.4 /jin radius. 

5. The value of C (x 0 /8, R/8) was converted to C (xo/6,0). Using the reference value 
(R/8) ief = 0 simplifies the comparison with the calculated values of C. 

Figure 13 plots the measured magnitude and phase of C (xo/6,0). The abscissa is the 
reduced displacement amplitude Xq/8 = r t i P 0o/^, where r t i p is distance of the oscillator’s 
tip from its torsion axis. The angular oscillations were sufficiently small {6 o < 0.05 rad) 
that the tip motion was always effectively rectilinear. Also shown are curves that represent 
the magnitude and phase of the calculated values of C {xq/8, 0) shown in Figure 11. The 
similarity of the measured and the calculated values is remarkable because the oscillator 
wire’s actual cross section is far from circular. Clearly, the condition 8 > R made the 
viscous force insensitive to the shape of the wire’s cross section. 

5 Concluding remarks 

Hydrodynamic scaling, numerical calculations, and experiment all show that the reduced 
force C on an oscillating cylinder depends chiefly on the scaling variable x$/8 when 8 > R. 
The dependences of C on the shape and size of the cylinder’s cross section are weak. Further 
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work building on these results might answer the following questions. 


• The scaling argument suggests that the force in a periodic flow held would scale as Xq/8 
even if the held were three-dimensional. Would this be supported by measurements or 
three-dimensional numerical calculations? 

• Eqs.(18) and (20) describe the dependence of C on R/6 over a limited range of x 0 /6. 
Over what ranges of R/6 and Xq/6 are these simple empirical forms valid? 

• What is C for a linear viscoelastic fluid? 

• Figure 1 shows that guitar strings and some insects oscillate near the boundary for 
two-dimensional how approximately defined by Xq/6 = 7. Why? 
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7 Appendix A: The moving reference frame approach 


This approach allows the hydrodynamic force on the cylinder to be calculated within a (non- 
inertial) reference frame that moves with the oscillating cylinder. In the laboratory (inertial) 
reference frame, the cylinder’s velocity is 


uq ( t ) = Re 


dt Xoe 


iujt 


—XqLU sin (cut) 


in the x-direction. At the cylinder surface, the fluid velocity is 


(29) 


u„ = (u,v) = (« 0 , 0) , (30) 

and far from the cylinder it is zero. Applying the velocity transformation u = u' + Uo 
creates a moving reference frame in which the cylinder is motionless and the fluid velocity 
far from the cylinder is — u 0 . Applying the same velocity transformation to the Navier-Stokes 
equation yields 


fin' 

P ( nt - u' ■ Vu' ) = -V (P+Po) + fiW, 


(31) 


where the time-dependent pressure is 


du 0 / 2 i i \ : 

po = P~Q^ X = x 0 cos (wt) X . 


(32) 
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Integrating p 0 over the surface of the cylinder relates the force per unit length F in the 
inertial reference frame to the force per unit length F' calculated in the non-inertial frame. 


F ( t ) = F' (t) — 7rR 2 pcu 2 x 0 cos (cut) 


( 33 ) 


Note that the force correction is proportional to tt R 2 p. More generally, one can show that the 
force correction for an oscillating body of arbitrary shape and fixed orientation is proportional 
to the fluid mass displaced by the body. 
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8 Appendix B: Inferred force on the oscillator’s tip 


This appendix describes how we inferred the force on the fastest moving wire, located at the 
oscillator’s tip, from the response of the entire oscillator. The oscillator response at large 
angular amplitude 9 o can be written in a form similar to that used at small amplitudes. The 
following is Eq.(15) in Reference [13] but generalized to large 9q. 





(34) 


(For clarity, Eq.(34) does not show the frequency-dependent corrections for anelasticity and 
the experimental electronics.) Here k and / 0 are the oscillator’s torsion spring constant and 
vacuum resonance frequency, respectively. The oscillator’s effective density is p s = pj ( irR 2 ), 
where p l is the wire’s measured linear density. The definition of p s assumes that the wire’s 
cross section is circular, unlike the actual wire. (See Figure 5.) Choosing R = 13.4 pm made 
the values of B osc measured at small amplitude similar to those calculated for a circular 
cylinder of the same radius. However, the size and shape of the assumed cross section is 
unimportant because it shifts only the magnitude of B osc . The shift does not affect the value 
of C because it is independent of oscillation amplitude. 

Inverting Eq.(34) yields B OSC (R/6,9 0 ), the dimensionless force-to-displacement ratio of 
the entire oscillator. Nonlinear hydrodynamics complicates the relation between 9 0 and 
Bosc because the experimental torque was a geometrically weighted sum of the local forces 
on individual oscillator wires. For large 9 o the local force on an individual wire was an 
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unknown function of xq, because the dependence of the amplitude function C on Xq/8 was 
unknown. 

Because the oscillator screen comprised some 200 squares, the screen was modeled as a 
continuous distribution of wires. Also, d 0 was sufficiently small (< 0.05 rad) that the local 
displacement Xq was proportional to 9 q. We thus approximated B osc (. R/6 , 9 o) by the integral 


B 


osc 



rs_/ 


3 

T* tip 






(35) 


Here, r measures distance from the torsion axis at r = 0, and the oscillator’s tip is located 
at r = r t i p . The product BC is the reduced force on a single screen wire that is located 
at r and oscillates with a linear displacement amplitude Xo = 9or. In the limit of small 
amplitude, (7 = 1, and Eq.(35) yields B osc (R/6, 0) = B (R/6), the function used in Eq.(6). 

To obtain the function that describes the oscillator’s tip we first take the derivative of 
Eq.(35) with respect to 9 0 . 



OSC 




d ( 9 0 r RY 
d9 0 yd’s) 

JL—c ( ®° r R ^\ 

9 0 dr V 8 ’ <5 ) 



(36a) 

(36b) 


Rearranging and integrating by parts yields 


B 



f xo R\ 

\8 ’ 8 ) 



9_o 

3 



( 37 ) 
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Eq.(37) gives the hydrodynamics of the oscillator’s tip in terms of the oscillator’s force/displacement 
ratio and its derivative with respect (){,. 
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9 Appendix C: Selected notation 


Table 6. Notation used in several places. 


A, A* 
ai,h 
B ( R/6 ) 

B cy i (R/6) 
B osc (f,d 0 ) 

C (xq/6, R/6) 
C 0 (xq/6) 

C<t> (xq/6) 

F 

-^small 

g (/, e 0 ) 

h, w 
KC 
L 

Nip 

R 

Uo 
U, V 

bin; bout 
X 0 

Nip (t) 

x,y 


Figures 1 and 2, two-dimensional flow regimes [10] 

Eq.(15), Fourier components of numerically calculated force 
Eq.(6), dimensionless force on arbitrary cylinder at small Xq 
Eq.(7), dimensionless force on circular cylinder at small Xo 
Eq.(34), dimensionless force on entire oscillator at amplitude 9 o 
Eq.(8), amplitude function 

Eq.(19), function that governs dependence of \C\ on Xq/8 
Eq.(21), function that governs dependence of <f) on xq/8 
Eq.(6), force per unit length on the cylinder or wire 
Eq.(6), force per unit length in the limit of small Xq 
Eqs.(26) and (34), experimental oscillator’s response 
Figure 7, height and half- width of computational domain 
Eq.(l), Keulegan-Carpenter number 
characteristic length for oscillating flow 

Eq.(25), distance between the tip and the torsion axis of the oscillator 

radius of the cylinder or characteristic radius of the wire 

velocity amplitude of the cylinder 

fluid velocity components in the x and y directions 

experimental input and output voltages 

displacement amplitude of the cylinder 

displacement of the wire at the oscillator tip 

coordinates in stationary reference frame 


P 

6 

V 

0o 

P 

P 


Eq.(2), Stokes viscous parameter 
Eq.(3), viscous penetration length 
viscosity 

angular displacement amplitude of the experimental oscillator 
density 

Eq.(17), phase of C (xq/8,R/8) 


Table 7. Modifications of some variables. 


modification 

example 

dimensionless 

velocity u* 

moving reference frame 

velocity u' 

Fourier transform 

force F 
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300 |— 



P 


Figure 1 . — Values of Keulegan-Carpenter number KC and viscous parameter 
p used in the calculations and measurements. The values for several familiar 
oscillators are also shown. The curve xo/5 = 7 is an extrapolation to small p 
of the two-dimensional flow regimes found by Tatsuno and Bearman’s flow 
visualization study [10], 



P 


Figure 2. — Data from Figure 3 of Tatsuno and Bearman’s flow visualization 
study [10], The curve xq/8 = 7 approximately bounds the two-dimensional 
flow regimes called A and A*. 
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Magnitude of C (xq/8, Rid) -n (R/S) 3/2 |B| 



igure 3. — Magnitude and phase of the hydrodynamic function B (Rid) for the present oscillator measured 
at small amplitudes [12], The analytical calculation is Stokes’ function B cy | (Rid). Choosing R = 13.4 pm 
for the oscillator’s effective radius brought the measured values near the calculated values. The similarity 
between the measured and calculated values shows that the oscillator’s hydrodynamics resembled that of 
an assembly of isolated circular cylinders, (upper) Magnitude scaled to reveal departures from the dominant 
(Rl 8) 3/2 behavior. The offset of the measured values is unimportant because it corresponds to an error of 
only 6% in the linear mass density of a single wire. (The nominal density was estimated by weighing a piece 
of screen different from that used to make the oscillator.) 



Figure 4. — Numerically calculated magnitude and phase of the amplitude function C (x^ld, Rid) as a function 
of reduced amplitude. The dependence of C on xo/8 is approximately independent of frequency even though 
R/d varies by a factor of V24. 
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Mag (C S h a p e )/mag (C c j rc | e ) 


0o 



Figure 5. — Schematic diagram of the oscillator (not to scale). 

The oscillator’s torsional displacement 0o caused the local 
displacement amplitude xo = rQ o, where r is the distance from 
the torsion axis. The inset shows the cross section of a typical 
oscillator wire superposed over the circular radius R = 13.4 pm 
used by the calculation. 
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Figure 6. — Numerical calculations at 1 Hz using three alternate cross sections showed that the shape of the cross 
section had only a small effect on the magnitude and phase of the amplitude function C (xo/8, Rl 8). The three 
shapes (shown in the inset) were a rectangle (8x30 pm), a rectangle with half inch ends (30 pm width, 4 pm radii), 
and a rectangle with quarter-circle ends (30 pm width, 8 pm radii). The results are normalized by those for the 
circle (1 3.4 pm radius) that was used in the main numerical calculations. The third shape resembled that of the 
experimental oscillator, and its results were used as a small correction of the experimental data. 
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reference frame 

Figure 7. — (a) Schematic illustration of the simplified flow problem of a circular cylinder oscillating in an incom- 
pressible incompressible fluid. In the laboratory (inertial) reference frame xOy the velocity of the cylinder is 
-xocosin(cof) in the x-direction. In the moving (non-inertial) reference frame x'Oy attached to the center of the 
oscillating cylinder the cylinder is motionless, (b) The domain and boundary conditions used in the two-dimensional 
numerical simulations. Use of the reference frame x'Oy made the cylinder stationary while the surrounding fluid 
oscillated sinusoidally in the x-direction. The lengths wand h defined the size of the computational domain. The 
domain included only y > 0 due to the assummed symmetry with respect to the x-axis. 



Figure 8. — Finite element mesh used in the validation test with 576 quadrilateral elements and 
a grading ratio q = 0.65. (a) Global view of the mesh bounded by -2000R < x < 2000 R and 
y < 2000R; (b) Enlarged view of the mesh grading near the cylinder in a local area bounded 
by -2 R < x < 2 R and y < 2 R. 


NASA/CR— 2007-2 1 5050 


34 




Magnitude 



r 


xq = 800 pm 



Figure 9. — Selected contours of the shear component of the rate-of-deformation tensor D xy indicated the concen- 
trations of viscous stress near the cylinder at 5 Hz. (Unit’s are s _1 .) Each map was calculated during the last 
cycle when cylinder was moving from right to left at maximum speed. Asymmetry is clearly visible at xo = 800 pm 
(x 0 /8 = 14.7). 




Figure 10. — The derivation of dC/dR was estimated from the difference C [xq/8, (13.4 pm)/8]-C [xq/8, (6.7 pm)/8]. 
(left) The normalized derivative of magnitude |C| with respect to Rl 8 was approximately described by a 
simple linear function of |C|. (right) The derivative of phase <|) with respect to R/S was approximately described 
by (|). (The curves are merely guides to the eye). 
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Figure 11. — Numerically calculated magnitude and phase of the amplitude function C (xq/8, 0) as a function of 
reduced amplitude. Eq. (22) was used to adjust the data to Rid = 0 which is equivalent to an extrapolation to 
0 Hz at constant xq/8. The curves are empirical description given by Eq. (23). 



Figure 12. — A single oscillator wire and its maximum displacement xo 
relative to the viscous penetration length 8. The values are for 3 Hz, 
which gave the largest value of xq/8. 




Xq/8 Xq/8 

Figure 13. — Measured magnitude and phase of the amplitude function C (xo/8, 0) as a function of reduced 
amplitude. The curves describe the calculated values derived from the data shown in figure 11. 
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